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Abstract 

We study the consequences of time variations in the scale of grand unification, My, when the 
Planck scale and the value of the unified coupling at the Planck scale are held fixed. We show that 
the relation between the variations of the low energy gauge couplings is highly model dependent. 
It is even possible, in principle, that the electromagnetic coupling a varies, but the strong cou- 
pling az does not (to leading approximation) . We investigate whether the interpretation of recent 
observations of quasar absorption lines in terms of time variation in a can be accounted for by 
time variation in M\]. Our formalism can be applied to any scenario where a time variation in an 
intermediate scale induces, through threshold corrections, time variations in the effective low scale 
couplings. 

PACS numbers: 12.10.Dni, 12.10.Kt, 98.80.Cq 
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I. INTRODUCTION 



Recent observations of quasar absorption lines have been interpreted as indicating vari- 
ation in the fine structure constant on cosmological time scales 0, 

— = (-0.57 ± 0.10) X 10"^ at 0.2 - 3.7. (1) 

a 

In the context of grand unified theories (GUTs), such a time variation in a may be related 
to time variations in other gauge couplings and, in particular, to variations in the strong 
coupling constant or, equivalently, in the QCD scale Aqcd- Such related variations affect 
various observables in interesting ways H, |, |, |, 0, |, |, ^ . 



Most previous works have assumed that the variation in a is induced by a variation in 
the unified coupling at the Planck scale, au{Mpi) H, |, |^, 0. The motivation for such 
studies comes mainly from string theories, where the value of au(^Pi) is determined by the 
value of the dilaton and other fields. We are here interested in another possibility, namely 
that the variation in a is induced by a variation in the GUT scale Mu, with Mpi and au(Afpi) 



held fixed. In ref. [|TT|, it is argued that if grand unification (in the sense of unification of 
the gauge interactions in a four- dimensional gauge group) arises within string theory, the 
GUT scale is likely to be a modulus.^ Such a situation can also arise in the framework of 
field theories, where symmetry breaking scales are determined by vacuum expectation values 
(VEVs) of scalar fields, but coupling constants are not. In either case, understanding why 
the modulus has the requisite properties is a difficult problem, which we will not address. 
We will show that the relation between 6a /a and 5Aqcd/Aqcd is highly model dependent. 
In particular, in contrast to the case that both are induced by 6a\j/a\j ^ 0, it may happen 
that 5AqcdMqcd ^ (though this is not the generic case). 

We are aware of two previous works that are closely related to ours. In ref. 0, the 
consequences of variations of Mu have been analyzed. It was assumed, however, that au(^u) 
remains fixed. It is then necessary that au(Mpi) also varies in a correlated way. We think 
that such a scenario is less plausible than the one that we investigate. In ref. [§, the 
consequences for the low energy couplings of variations in the mass of particles have been 
analyzed. Conceptually, our work takes a similar direction. In their discussion of grand 
unified theories, ref. [Bl examine a set of models in which only the masses of a subset of 



^ Of course, in string theory, coupling unification does not necessarily require a unified gauge group. 
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GUT mass particles vary in time; in particular, only the masses of color-neutral, charged 
particles shift. In contrast, our basic assumption is that all heavy masses (of GUT gauge 
bosons, fermions and scalars) are proportional to a single breaking scale, Mu, and it is the 
variation of this scale that we study.^ 

In this work we focus on the resulting relations between the low energy couplings. We 
do not concern ourselves here with the origin of the scale variation, nor with the connection 



to higher scale physics. Discussions of these topics can be found in g [T^. In addition, 
there are fundamental issues we will also not address. As pointed out in it is difficult 
to understand, in the framework of local quantum field theory, how couplings could vary by 
so large an amount, without an enormous variation in the vacuum energy. This problem is 
not significantly different if it is changes in the unification scale which are responsible for 
this variation than if the coupling at the Planck or unification scale changes. 

The plan of this paper is as follows. In section || we present the formalism for the variation 
of low energy couplings induced by variations in the scale of threshold corrections. In 



section |IIT| we obtain the relations between the variation of relevant cosmological observables 
and argue that, within our framework, they depend on a single parameter. In section ^ 
we apply the formalism to the experimental data. We use it to constrain the parameter of 
section In section |V|we test specific models of grand unified theories with this constraint. 
We show that a variation in the GUT scale is unlikely to explain the claimed variation in 
6a /a. In section |VT|we explain how our formalism and results apply much more generically, 
to any model where there is a time variation in a scale where threshold corrections take 
place. We also mention some further subleading effects which might be non-negligible. We 



conclude in section VII. 



II. FORMALISM 



We assume that the variation in couplings is due to a variation in a single intermediate 
scale, Mu, where GUT threshold corrections take place. The gauge couplings ai{Q) at a 



^ To understand the difFerenee, one may think of the consequences of a variation in the VEV of the single 
Higgs doublet of the Standard Model, which would shift all masses, or a shift in the VEV of a hypothetical 
Higgs triplet, which would shift only neutrino masses. 
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«-'(Q)=«u'(Mu) + ^ln(^- 
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(2) 



scale Q < Mu are given to one loop by the RGE solution 

Q 

where au(Mu) is the unified coupling at the scale M\j and hi is the beta function coefficient 
(between Mu and Q) for Oj. The scale Mu is related to the mass scale of heavy degrees 
of freedom that are integrated out, particularly the heavy GUT gauge bosons. Since the 
particle content of the theory below and above Mu is different, the RGE coefficient hi changes 
at Mu- In particular, one obtains 



«u^(Mu) =«u^(Mpi) + ^ln ( 

/vr \ 



where h\j is the beta function coefficient for the unified group. Combining (^) and we 
write 

Now we allow the GUT scale, Mu, to change by an amount of 5Mu, while holding au(Mpi) 
and Mpi fixed. The resulting variation in the low scale couplings can be calculated from 
(a similar derivation, for oem, appears in [§): 

5ai{Q) Ahi6My 
ai[Q) Z-K Mu 

where 

Ahi = hu-hi. (6) 

Using d^), we can now relate the variation in low scale parameters to the variation 6 My. 
We focus on two parameters, the electromagnetic coupling a and the QCD scale Aqcd^ 

a-' = ^a^\0) + a^\0) , (7) 

A«co^Mr'™r™-exp(-^). (8, 

Neglecting effects of threshold corrections below Mu (which are discussed later), we obtain 
the relations: 

Aqcd 9 Mu ■ ^ ^ 
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Finally, we use @ and (0) to relate the variation of the two parameters: 

i^Aqcd _ 27r A63 5<y_ 
Aqcd ~ 9a IA61 + A62 a 

^ i? ^. (11) 

a 

The factor R ^ 5 Ab''+Ab P^^Y^ ^ crucial role in our analysis. As we show later, its value is 
highly GUT- model dependent. 



III. OBSERVABLES 



In this section, we use (|lID to obtain expressions for the variation of several cosmological 
observables. The primordial ^He abundance is given by [1^, |15 



_ 2(n„/np)Ns 
1 + (n„/np)Ns 

where the ratio of the number density of the neutron to that of the proton at nucleosynthesis 
is given by 

(nn/np)Ns ^ 0.8 x e'^/^^ . (13) 

Here ~ 0.8 MeV is the decoupling temperature and Q is the mass difference between 
the neutron and the proton. Since To is a function of Mpi and Gp, it is independent of the 
Standard Model gauge couplings and we therefore neglect its variation. The mass difference 



Q, on the other hand, depends on both a and Aqcd ||T6 



Q - 2.05 MeV + C«Aqcd , (14) 

where C is a dimensionless order one parameter. Assuming variations in a and Aqcd, the 
value of Q is shifted [eq. from its present-day value of Q — 1.29 MeV, and consequently 
i^n/ rip) ]<fs would be shifted [eq. ([T3|) ] from its 'standard' value of (n„/?T,p)Ns — 7- We thus 
obtain the change in F4 compared to its value calculated with the present values of the 
coupling constants: 



SYj _ 1 2.05 MeV - Q /6a ^ 6Aqcd\ 

Y4 l + (n„/np)Ns Td \a Aqcd/ 



0.8(1 + m — . 

a 



(15) 
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Two other relevant observables are X and Y. The X parameter gives the ratio between 
the hyperfine 21cm neutral hydrogen absorption transition to an optical resonance transi- 
tion fl^: 

A-.„^,,(!^), (16) 

where rUe is the electron mass, rup is the proton mass, and Qp is the proton gyromagnetic 
moment. The Y parameter determines the molecular hydrogen transition in the early uni- 
verse fO, IT 



F = ^ . (17) 



The proton mass is proportional, at first order, to Aqcd 



6mp ^ 5Aqcd _ j^Sa ^^^^ 



rrip Aqcd a 
Assuming that the variations of rrie and Qp are negligible, we obtain: 



^^2^-^^(2-i?)^, (19) 
X a m„ a 



-^R-. (20) 



IV. CONSTRAINING R WITH EXPERIMENTAL DATA 

Values for the observed variation of various observables were obtained through cosmolog- 
ical and nuclear studies. Table | summarizes these values (the details can be found in the 
quoted references) . Using the formalism of the previous section, we relate the allowed range 
for each observable in Table | to an (in general, -R-dependent) allowed variation in a. The 
results are also given in Table |. 

For a given value of R, one obtains from Table | the implied variation of a at different 
z^s. The most interesting result, however, arises from requiring consistency of the three 
observables related to the intermediate redshift values. Note that 6Y/Y and 6X/X are 
consistent with each other for any value of R, because, at the 2a level, they allow 6a/a = 0, 
for which the value of R becomes irrelevant. The constraints become nontrivial when one 
considers the claimed signal of 6a/ a [eq. (|l])]. The requirement of consistency between 
the observed 6Y/Y at z ^ 2.3-3 and the observed 6a/a at z ~ 0.2-3.7 favors a range 
in R that is mostly negative in sign and not very large in magnitude, R = O{—10). The 
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TABLE I: Cosmological and nuclear values for the variation of various observables and their im- 
plications for 5a/ a. 



Constraints 



Redshift z 



Reference 



Implied constraint on 5a/ a 



^ = (-6±7)xl0-3 

— = (-0.57 ±0.10) X 10" 

a 



5Y_ 

'V 

5_X_ 

— = (-0.36 ± 1.44) X 10" 

a 



(5.7 ±3.8) X 10"^ 
(0.7 ± 1.1) X 10~^ 



lOio 
0.2-3.7 
2.3-3 
1.8 
0.1 





SI 



(-7 ±9) X 10-3/(1 + ^) 
(-0.57 ±0.10) X 10-5 
(5.7 ±3.8) X lO^V-R 
(0.7 ± 1.1) X 10-V(2 - 
(-0.36 ± 1.44) X 10-*^ 



requirement of consistency between the observed bX/X at z ^ 1.8 and the observed 8a/ a 
at z ~ 0.2-3.7 favors a range in R that is mostly positive in sign and, again, not very large 
in magnitude, R = 0{+3). The most interesting results follows from fitting simultaneously 
all three observables. The result of a analysis is that the following range of R is favored: 



-1 < i? < +6 



(21) 



Values of R outside this range have a probability that is ^ 5% to yield the observed results. 
(Note that even the best fit point, with R ~ +2, has only a probability ~ 15% to yield the 
observed results.) 

This result can be understood qualitatively: When |_R| is large, a value for 6a/ a implies 
an even larger value for (5Aqcd/Aqcd- This larger variation dominates the shift of nucleon 
masses. Since Y is proportional to the proton mass while X is inversely propotional to it, 
their relative variations must be of opposite signs, 6Y6X < 0. The data, however, indicate 
that their variations are both biased to be positive. Consequently, only one of them can 
be consistent with a non-zero variation of a. (Which one depends on the sign of R.) This 
situation is avoided only in case that the variation in Aqcd does not dominate 6X/ X or, in 
other words, in case of a small R. 

To summarize, eq. (PTT) gives the consistent range for the parameter R in theories where 
the variation of the low energy coupling constants is dominated by threshold effects at a 
single varying scale. 
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V. IMPLICATIONS FOR VARIOUS GUT MODELS 



In this section we apply our results to specific models. We focus on the minimal super- 
symmetric standard model (MSSM) embedded in a grand unified theory (GUT), in which 
the varying scale is the scale of breaking of the unified group, M\j. In general, this is the 
mass scale of the GUT gauge supermultiplets, and we assume that it characterizes also the 
masses of all the heavy chiral supermultiplets. The one loop beta function coefficients below 
Mu are those of the MSSM: 

bi = -, h = l, h = -3. (22) 

5 

Above the GUT scale, there is a single coefficient common to all gauge couplings, by. Thus, 
the ratio R of eq. (|ll]) is given by 

27T bv + 3 



R=^ .7 . (23) 
9a |6u-12 ^ ' 



Before we discuss specific GUT models, let us make some general observations: 

1. For asymptotically large (positive or negative) values of b\j, we obtain 

|&u|>l =^ i?^^~+36. (24) 

Large negative values of bu are not relevant, as demonstrated by the very conservative lower 
bounds set by examining the contributions from the gauge supermultiplet and from the 
minimal matter representations that are necessary to accommodate the MSSM fields: 

bij[SU{h)] > -8, by[SO{10)] > -17. (25) 

On the other hand, realistic models often have a large, positive b\j and predict R = (9(40). 

2. Negative values of R are achieved only for a very limited range of b\j: 

< for - 3 < 6u < (26) 



3. For a hypothetical value of b\j = ^, a does not vary (to the approximation that we 
use): 



bv = l ^ |i?|^oo. (27) 

This result explains the change of sign of R between models with bu = 4 and bu = 5, and 
the very large magnitude of R in these models, R — 502(+574) for bu = 4(5). 
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4. For bu = bs = —3, Aqcd does not vary (to the approximation that we use): 

6u = -3 =^ R = 0. (28) 

In general, the (9(100) factor in eq. (pAf), which is due to the exponential dependence of 
Aqcd on as, leads to a large R. In order to overcome this factor and obtain \R\ ~ 0{1), the 
threshold correction of 0:3 needs to be highly suppressed, i.e., A63 < O.OIA61 2. For MSSM 
GUT models with varying Mqut, this would imply —3.5 ^ &u ~ —2.5. In particular, when 
bu is shifted by one unit from the special value of —3, one already obtains rather large values 
of |-R|, R = +4.2(— 5.5) for by = — 4(— 2). We learn that, if R is to be within the range of 
eq. (pID , the only acceptable integer values for b\j are —3 and —4. 

We now move on to discuss several examples of specific GUT models. The bu coefficients 
and the resulting i?-factors of various models are given in Table |I[ Each model is defined 
by its representation content. (Note that, for example, gives the number of SU(5) funda- 
mentals plus the number of antifundamentals.) All SU(5) models have the quark and lepton 
fields in three generations of 10 + 5 and the Higgs fields related to electroweak symmetry 
breaking in 5 + 5. The minimal model has, in addition, a single 24 to break SU{5) — * Gsm- 
This model is, however, excluded by combining constraints from coupling unification and 
from proton decay (It may still be viable with a special fiavor structure of the super- 
symmetric mixing matrices [^.) A viable, though fine-tuned, model can be constructed by 
adding another 5 + 5 pair To naturally induce doublet-triplet splitting, the 24 has 



to be replaced by, at least, 50 + 50 + 75 |2|, |3> Hf. 

All SO(IO) models have the quark and lepton fields in three generations of 16. The min- 
imal Higgs sector has 45 + 16 + 16 to break 5*0(10) — > Gsm and a single 10 for electroweak 



symmetry breaking p9[. Note that we assume here that the breaking of 5*0(10) Gsm is 
done in one step, and that the VEVs of all GUT-breaking Higgs fields vary together. The 
investigation of a two-step breaking, 5*0(10) SU{5) Gsm, with independent variation 
of each scale, will be presented in future work.'^ Models where doublet-triplet splitting is 
achieved naturally require an additional 10 and either an additional 16 + 16 pair or a 



^ In case that the scale of S'O(IO) — > SU{5) breaking varies but the scale of SU{5) Gsm breaking does 
not, the result is effectively a variation in the unified coupling constant at the SU{5) breaking scale. This 
result shows how our formalism can be applied also to models of time variation in ay: we should take 
Abi = A62 ~ A63, which gives i? ~ 36. 
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more complicated Higgs sector that includes several 45 and 54 multiplets In each of 
these models, one could have a 126 + 126 pair in place of a 16 + 16 pair. Such a replacement 
increases by by 66 and gives R = 0(40). 



TABLE II: Various GUT 


models, their beta funct 


ion coefficient bu and the resulting 


factor R. 


SU(5) 








Tic;n Tire; 


feu 


R 




5 


3 


1 




-3 







7 


3 


1 




-2 


-5.5 




5 


3 




2 1 


+52 


+41.5 


SO(IO) 




ni6 


"■45 


7154 


feu 


R 




1 


5 


1 




-5 


+7.6 




2 


7 


1 







-23.9 




2 


5 


3 


2 


+36 


+44.4 



Other models of interest are reviewed in, for example, ref. o2 and lead to similar results. 
We conclude that most GUT models give ^ 1. A remarkable exception is given by 
models with by = —3. This possibility is realized in the minimal SU{5) model which, 
without a very special flavor structure, is, however, phenomenologically excluded. Other 
models with by = —3 can be constructed within the framework of larger groups. There is, 
however, little motivation to consider such a scenario (which would require fine-tuning for 
the doublet-triplet splitting) for the sole purpose of deriving the desired by. We conclude 
that, within the framework of realistic GUT models that avoid fine-tuning, our mechanism 
predicts \R\ ^ 1. 

When this result is confronted with eq. (pi]), one is led to the conclusion that supersym- 
metric GUT models with a varying GUT scale are unlikely candidates to explain a variation 
in a of the size given in eq. (|1|) . 



VI. GENERALIZATION AND LIMITATIONS 



Our formalism and many of our results have a much broader applicapibility than the My- 
variation in GUT models that we have focussed on. The formalism applies to any variation 
in a single intermediate scale where threshold corrections take place, fith- Suppose that the 
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couplings are held fixed at some fixed high energy scale /io- Define bi to be the one- loop 
beta function coefficients between fith and a low scale Q, while b'^ are the corresponding 
coefficients between /xq and /Xth- We now define, as a generalization of eq. (||), 

Ab, = b[- h. (29) 

With this definition, eq. (|^) for 5a/ a, eq. (|TUp for (5Aqcd/Aqcd, and eq. (|TT]) for the ratio 
i?, still apply.^ Similarly, the constraints on 5a /a in Table H and the allowed range for R in 
eq. (|2l|) hold in these more general circumstances. 



An important comment is now in order: In principle, the variation (5/ith may induce a 
variation in all couplings at lower scales. In the most general case, this includes a variation 
in the masses of all the particles which are lighter than /ith. A variation of these masses 
could be induced by two sources: 

1. The parameters that determine a particle mass, such as Yukawa couplings, may have 
threshold corrections at /ith- In this variation analogous to (^) results. To 
calculate this effect, one should use the (linearized) beta function of these parameters 
below and above the scale /ith • 

2. The variation in the gauge couplings, below /ith, leads to variations in the low scale 
masses. To calculate this effect, one should use the RGE of the mass parameters below 

/^th- 

When running the gauge couplings from /ith to some low scale Q, all particles with masses 
Q < m < /ith are integrated out. For example, in an MSSM GUT model, with Q ~ rrib 
and /ith = ^GUT, we integrate out all supersymmetric particles at a scale Msusy as well 
as heavy SM particles (t, Z, H^^). This integration out would usually introduce additional 
threshold corrections, which have small effects on the parameters at rrih. 

We are mainly interested in small variations of a and Aqcd- As explained above, the 
variation 5/ith induces variations in all intermediate scales below /ith- But these variations 
induce, by threshold corrections at these scales, further variation in gauge couplings. The 



^ One has to be careful about the normahzation of ai in non-GUT models. For example, if we take the 
conventional definition of U{1)y in the Standard Model, the factors of I in eqs. (|), (|) and have to 
be omitted, while for eq. (Ea) one should use bi = 11. 
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whole process may be pictured as a 'chain reaction', with a final effect that could be signif- 
icant, even if the individual variations at each intermediate scale are small. We postpone 
the detailed study of this effect to future work. 

Here, we would only like to emphasize that when R is large, as is the case in most of 
the GUT models that we have investigated, the variation of Aqcd is the dominant effect at 
low energies, and it is well justified to ignore variations in couplings other than the gauge 
couplings. For small R, however, the modification due to threshold corrections of other 
parameters can become important. Thus, when we think of i? = models, one should not 
conclude that a would vary with Aqcd remaining strictly fixed. Very likely, Aqcd varies 
due to the effects of lower thresholds and of higher loops. But the effect is expected to be 
small [R = 0{1)] and our qualitative conclusions would not change. 



VII. CONCLUSIONS 



A variation in a physical scale in which threshold corrections take place leads to varia- 
tions in low scale observables. Current experimental data constrain the relations between 
the threshold corrections of the three different gauge couplings. In particular, consistency 
between the claimed variation 6a /a [eq. (0)] and the allowed ranges for 6X/X and 6Y/Y 
[Table H give 

. 27r A63 . 

where A6j is the difference of the one-loop beta function coefficient above and below the 
varying scale. In other words, the difference between the beta function coefficients of the 
strong coupling 0:3 above and below the varying scale, should be suppressed by a factor of 
(9(0.01) compared with those of the other two couplings, ai and 0^2. 

We focussed our investigation on the framework where the MSSM is embedded in GUT 
models and the varying scale is Mqut- We demonstrated that the relation between the 
variation of Aqcd and that of a, that is, the ratio R = is highly model dependent. 

For example, we can construct SU(5) models with R as high as (9(+600) and as low as 
500), and a model where R = 0. (The latter is the minimal SU(5) model, and the 
former have, respectively, eight and seven 5 + 5 pairs added to the minimal model.) We 
argued, however, that it is difficult to obtain consistency with the requirement (|30|) in realistic 
supersymmetric GUT models in our framework. 
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